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The stability and properties of the monovacancy and the divacancy in single-walled carbon nan- 
otubes (CNTs) are addressed by spin-density functional calculations. We study these defects in 
four nanotubes, the armchair (6,6) and (8,8) and the zigzag (10,0) and (14,0), which have diameters 
of about 8 and 11 A, respectively. We also study different defect concentrations along the tube 
axis by increasing the supercell in this direction in order to have one defect every 13 and 26 A of 
CNT length. Our results show that in the equilibrium geometry CNTs with a monovacancy exhibit 
ferromagnetism with magnetic moments ranging from 0.3 to 0.8 fiB- On the other hand, CNTs with 
a divacancy do not exhibit magnetism due to the full reconstruction around the defect where all 
C atoms are three coordinated. We observe that the presence of a monovacancy does not change 
drastically the CNT electronic properties, preserving their corresponding metallic or semiconduct- 
ing character. However, for the divacancy both armchair and zigzag CNT become semiconductors 
exhibiting a energy gap of about 0.15 eV. 



I. INTRODUCTION 



The study of radiation-induced defects in car- 
bon nanostructures have become considerably impor- 
tant since the recent experimental results of room- 
temperature magnetic ordering in proton-irradiated 
graphite, which exclude metal impuritiesi, as well as 
the observation of ferromagnetism in pressure-induced 
polymerized fullerenes 2,3 . Magnetic ordering in de- 
fective graphite and fullerenes have been also pre- 
dicted theoretically 4,5 , suggesting the possibility to re- 
alize nanoscale magnets with potential applications in 
spintronics and biocompatible magnetic materials. On 
the other hand, magnetism in single-walled carbon 
nanotubes (CNTs) appears to be more complex than 
graphite and fullerenes owing the curvature, chirality and 
different electronic characteristics of the tubes. Single- 
walled CNTs, which are typically 10 A in diameter, can 
be metallic or semiconducting depending on the way that 
a graphene sheet is rolled up, which is characterized by 
the chiral indices (n, m). Metallic tubes occur if n — m is 
divisible by 3, otherwise the tubes are semiconducting 6 . 
Tubes with indices (n, 0) and (n, n) are termed zigzag 
and armchair, respectively, which is related to the ar- 
rangement of carbon atoms around the tube. 

Experimental knowledge concerning defect-induced 
magnetism in CNTs is still lacking. The difficult to mea- 
sure magnetic ordering in nanotube samples resides in 
the presence of metal catalysts which are necessary to 
produce single-walled tubes. Theoretical investigations 
on open ended zigzag nanotubes are found to exhibit 
energetically favorable ferromagnetic spin configurations 
which are associated to unsaturated dangling bonds at 
zigzag edges^. Single vacancies in small-diameters CNTs, 
ranging from 4 to 8 A, have been recently addressed by 
first-principles calculations 7 . Although relative high con- 
centrations of vacancy or small supercell sizes containing 
a vacancy have been considered in this work, it shows 
that only metallic nanotubes with a single vacancy in 



their ground state equilibrium geometries would exhibit 
ferromagnetism. 



II. THEORETICAL APPROACH 

The calculations were performed in the framework 
of the spin-polarized density functional theory 8 , using 
a basis set of strictly-localized numerical pseudoatomic 
orbitals as implemented in the SIESTA code 9 . The 
exchange-correlation energy is calculated within the gen- 
eralized gradient approximation's. Standard norm- 
conserving pseudopotentials 11 in their separable formic 
are used to describe the electron-ion interaction. We use 
a double-^ singly-polarized basis set. We study the va- 
cancies in both armchair and zigzag CNTs with different 
diameters. The armchair CNTs have chiral indices of 
(6,6) and (8,8), whereas the zigzag ones have indices of 
(10,0) and (14,0). For the (6,6) and (10,0) CNTs which 
have diameters of about 8 A, we use unit cells containing 
120 atoms. For the (8,8) and (14,0) CNTs with diameters 
of about 11 A, the supercells contain 160 and 168 atoms, 
respectively. For both 8 and 11 A diameter CNTs the 
supercell has a length of L « 13 A. This means that we 
are calculating one defect every 13 A of CNT length. For 
the Brillouin zone sampling we use six k points along the 
CNT axis 13 . To ensure negligible interaction between 
tubes we imposed a vacuum region of 12 A. With the 
above supercells, the distance between vacancy images 
along the tube axis is about 13 A. This means that we 
are simulating a CNT with a linear concentration of de- 
fects with one defect every 13 A of CNT length. A longer 
supercell with length of 2L w 26 A are also consider to 
study the effect a lower defect density in the CNT elec- 
tronic structure. To do that we use (6,6) and (10,0) nan- 
otubes in supercell containing up to 240 atoms and the T 
point for the BZ sampling. The positions of all atoms in 
the nanotubes were relaxed using the conjugated gradi- 
ent algorithm until the force components become smaller 
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FIG. 1: Equilibrium geometry of the monovacancy and 
the divacancy in (6,6) and (10,0) CNTs. (a) (6,6)+lV, (b) 
(10,0) + 1V, (c) (6,6)+2V, (d) (10,0)+2V. Painted balls dis- 
play the C atoms around the defect. 



than 0.05 eV/A. 



III. THE MONOVACANCY IN CARBON 
NANOTUBES 

Figures 1(a) and 1(b) show respectively the equilib- 
rium geometries of the (6,6) and (10,0) CNTs contain- 
ing the monovacancy (hereafter referred to as (6,6) +1V 
and (10,0)+1V). We find that the neighboring C atoms 
partially reconstruct around the defect, forming a pen- 
tagon bonding structure, leaving a C atom twofold co- 
ordinated, which is projected slightly off to the nan- 
otube surface. We also find met ast able configurations 
for (6,6)+lV and (10,0)+1V, 10.5 and 12.1 meV/atom 
higher in energy that their respective ground state (GS) 
geometries. In these configurations the pentagon struc- 
ture plus the twofold coordinated C atom are also found 
but rotated clockwise about 60° with respect to the equi- 
librium geometries. A third met ast able structure is found 
for (10,0)+1V, which is 17.2 meV/atom higher in energy 
than the GS geometry. Here, the C atoms neighboring to 
the vacancy do not form bond, remaining twofold coor- 
dinated. This geometry is known as 3db (three dangling 
bonds) structure. 



Similar results for the GS and metastable geometries 
are found for the vacancy in the 11 A diameter CNTs 
(8,8) and (14,0). The metastable geometries of (8,8)+lV 
and (14,0)+1V are 7.3 and 6.7 meV/atom higher in en- 
ergy that the respective GS geometries. The smaller en- 
ergy difference with respect to the 8 A diameter CNTs is 
attributed to a decrease in the strain energy due to the 
lower curvature of the (14,0) and (8,8) tubes. Thus, ex- 
trapolating to higher-diameter tubes, we can infer that 
the metastable geometries would tend to disappear. It 
is interesting to note that the 3db metastable structure 
of (10,0)+1V is not found in (14,0)+1V, suggesting that 
the 3db geometry exits due to the higher strain in the 
(10,0) tube, hindering further relaxations. 

We calculate the formation energy of a vacancy in our 
CNTs using, 

Ef rm(V) = E tot (CNT + V) - E tot (CNT) + \i c (1) 

where £^(CNT+V) and £ tot (CNT) are the total energy 
of the tube containing the vacancy and the perfect tube, 
respectively. \iq is the carbon chemical potential which 
is calculated as the total energy per atom in the perfect 
tube. Our results are listed in Table I where we also in- 
clude vacancy calculations performed in (6,6) and (10,0) 
CNTs considering a supercell with twice the length of the 
one previously used. We find that the vacancy in the 8 A 
diameter CNTs have similar formation energies, ranging 
from 5.65 to 5.85 eV, including the longer supercell. The 
metastable geometries also show the same tendency, with 
formation energies of about 7 eV. The above suggests a 
minor influence between images vacancies for supercell 
with lengths up to L ^13 A. For the vacancy in CNTs of 
11 A of diameter we find lower formation energies than 
those of the 8 A ones, where the defective tube (8,8)+lV 
[(14,0)+1V] has a formation energy of about 1 eV [0.5 eV] 



TABLE I: Formation energy (Ef orm ) and magnetic moment 
(m) for the monovacancy and the divacancy in armchair and 
zigzag CNTs, considering different supercell length (L) or va- 
cancy concentrations, dc-c is the equilibrium distance be- 
tween C atoms that approach each other to form the pentagon 
bonding structure. In parenthesis are shown our results for 
the metastable geometries. 



System 


L(A) 


m (ijlb) 


Ef orm 


(eV) 


dc-c 


:(A) 


(6,6)+lV 


12.46 


0.82 (1.1) 


5.75 


(7. 


,00) 


1.56 


(1.93) 


(6,6)+lV 


24.92 


0.63 


5.85 






1.57 




(10,0) + 1V 


12.96 


0.32 (1.0) 


5.67 


(7. 


.11) 


1.54 


(1.72) 


(10,0)+1V 


25.92 


0.49 


5.65 






1.53 




(8,8)+lV 


12.46 


0.77 (1.1) 


4.65 


(5. 


,81) 


1.58 


(1.98) 


(14,0)+1V 


12.96 


0.53 (1.1) 


5.20 


(6. 


.32) 


1.58 


(1.74) 


(6,6)+2V 


12.46 


0.0 (0.0) 


4.24 


(7. 


.84) 


1.53 


(1.78) 


(6,6)+2V 


24.92 


0.0 


4.17 






1.52 




(10,0)+2V 


12.96 


0.0 (0.0) 


3.90 


(6. 


.65) 


1.50 


(1.65) 


(10,0)+2V 


25.92 


0.0 


4.06 






1.49 




(8,8)+2V 


12.46 


0.0 


2.09 






1.53 




(14,0)+2V 


12.96 


0.0 


2.33 






1.50 
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FIG. 2: Spin-density isosurfaces for the monovacancy in 8 A 
diameter CNTs. (a) (6,6) + lV and (b) (10,0) + 1V. The iso- 
surfaces correspond to a charge of 0.015 e/A 3 . 



lower than that found for (6,6)+lV [(10,0)+1V]. The de- 
crease in the formation energy is attributed to the de- 
crease in the strain energy due to the small curvature of 
the larger tubes. 

Table I shows our results for the magnetic moment of 
the defective CNTs. We find that all nanotubes under 
consideration exhibit ferromagnetic ordering. The lower- 
energy structures of (6,6)+lV and (10,0)+1V have mag- 
netic moments of 0.82 and 0.32 /iB, respectively, whereas 
for the metastable configuration the magnetic moment 
are closer to 1 /iB- Clearly, the magnetism is due to the 
under-coordinated C atom which has a localized unpaired 
spin, as shown in Fig. 2. However, as a C atom with a 
dangling bond has a magnetic moment of 1 /iB , the lower 
magnetic moments found in (6,6)+lV and (10,0)+1V are 
due to a redistribution of charge around the defect. This 
can be checked by looking at the equilibrium distance 
between the C atoms that form the pentagon structure, 
shown in Table I. As longer this distance is, larger is 
the magnetic moment. This means that once a vacancy 
is created, two of the three dangling bonds containing 
a spin rehybridize forming the pentagon structure. The 
weakness of this bond depend on the curvature strain of 
the tube. An extreme case occurs in the 3db configu- 
ration of (10,0)+1V. Here no bond between the under- 
coordinated C atoms is formed, resulting in three dan- 
gling bonds with magnetic moment of 1.4 /iB- Consid- 
ering lower concentration of vacancies, that is doubling 
the supercell length, we find that magnetic moments on 
(6,6)+lV and (10,0)+1V are 0.63 and 0.49 /i B , respec- 
tively. Thus, the magnetic moment tend to stabilize ap- 
proximately around the average value 0.6 /iB for more 
diluted monovacancy concentrations. The same behavior 
is observed for the tubes with higher diameters where the 
magnetic moment on (8,8) + lV and (14,0)+1V are 0.77 
and 0.53 /iB, respectively. This stabilization appears to 
be related to a decrease in the strain energy on the tube 
surface. Thus, we may infer that the magnetization due 
to a single vacancy would be present in a wide variety of 
CNTs. 
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FIG. 3: Spin-resolved density of states (DOS) for the (6,6) 
CNT with a monovacancy. (a) The perfect (6,6) CNT (b) 
(6,6)+lV calculated with a supercell of 12.46 A length, (c) 
(6,6)+lV calculated with a supercell of 24.96 A length. In 
(a) and (b), upper and lower panels correspond to majority 
and minority spin configurations, respectively. Dashed lines 
indicate the position of the Fermi level (Ef). 



Figures 3(b) and 3(c) show the spin-resolved density 
of states for (6,6)+lV CNT calculated with both super- 
cells length (12.46 and 24.96 A respectively), which are 
compared with density of states of perfect (6,6) CNT 
[Fig. 3(a)]. We observe that the (6,6) armchair CNT 
preserves its metallic character for different vacancy con- 
centrations. We also note that the difference between 
majority (upper panel) and minority (lower panel) spin 
configurations in Fig. 3(b) and 3(c) clearly demonstrate 
the magnetism in the (6,6)+lV system. Similar results 
are found for (8,8)+lV CNT which suggests that these 
properties would be also found in larger diameter arm- 
chair CNTs. Figures 4(b) and 4(c) show the spin-resolved 
density of states for (10,0)+1V CNT calculated with su- 
percells of 12.46 and 24.96 A length, and the density of 
states of perfect (10.0) CNT [Fig. 4(a)]. As we can see the 
(10,0) semiconducting zigzag CNTs maintain this prop- 
erty when a C atoms is removed. For higher vacancy 
concentration the gap energy is about 0.1 eV mainly due 
to a state located just above the Fermi energy. For lower 
vacancy concentration this state is less dispersive increas- 
ing the gap energy up to 0.2 eV. The presence of this 
states give to the (10,0)+1V CNT an acceptor character. 
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FIG. 5: Density of states for the (6,6) CNT with a diva- 
cancy. (a) Perfect (6,6) CNT (b) (6,6)+2V calculated with a 
supercell of 12.46 A length, (c) (6,6)+2V calculated with a 
supercell of 24.96 A length. Dashed lines indicate the position 
of the Fermi Level (Ef). 



FIG. 4: Spin-resolved density of states for the (10,0) CNT 
with a monovacancy. (a) The perfect (10,0) CNT (b) 
(10,0) + 1V calculated with a supercell of 12.96 A length, 
(c) (10,0) + 1V calculated with a supercell of 25.92 A length. 
Dashed lines indicate the position of the Fermi Level (Ef). 



IV. THE DIVACANCY IN CARBON 
NANOTUBES 

Figures 1(c) and 1(d) display the GS geometries for the 
divacancy in (6,6) and (10,0) CNTs (hereafter referred 
to as (6,6)+2V and (10,0)+2V). We also find metastable 
geometries for (6,6)+2V and (10,0)+2V, which are 10.5 
and 23.3 meV/atom higher in energy than the respec- 
tive GS geometries. We observe that once a divacancy is 
created, the undercoordinated C atoms spontaneously re- 
construct around the defect forming an octagon with two 
adjacent pentagon in opposite directions. For (6,6)+2V 
the pentagons point along a line of about 30° with re- 
spect to the tube axis, whereas in the metastable struc- 
ture the pentagons are aligned normal to the tube axis. 
For (10,0)+2V the pentagons align along the tube axis 
for the GS geometry and along a line of about 60° with 
respect the tube axis. Similar reconstructed geometries 
for (8,8)+2V and (14,0)+2V CNTs where found, suggest- 
ing that these reconstructions should be found in larger 
diameters CNTs. For the 8 A (6,6) and (10,0) CNTs 
the divacancy formation energies are found to be 4.2 and 



4.0 eV, respectively. Using longer supercell we note a 
small change in the formation energies suggesting a neg- 
ligible interaction between the defect and their images 
in neighboring supercell. For the 11 A (8,8) and (14,0) 
CNTs the formation energies decrease to 2.1 and 2.3 eV, 
respectively, about half the value found in the 8 A CNTs. 
This suggests that the curvature strain becomes impor- 
tant for the divacancy formation in small diameter CNTs. 
In Table I we also show the bond length between the C 
atoms that form the pentagon structure. The smaller 
bond length, of about 1.5 A, are found for the divacancy 
in the zigzag (10,0) and (14,0) CNTs. As a consequence 
of the full reconstruction of the divacancy all C atoms 
around the octagon are threefold coordinated. Thus, no 
unpaired a electron is found, implying that carbon nan- 
otubes with this defect do not exhibit magnetism. There- 
fore, we can infer that the origin of the magnetism in 
nanotubes and graphite are mainly due to unpaired a 
electrons which are localized. On the other hand, the 
delocalized 7r electrons would have a negligible contribu- 
tion because once an unpaired a electron reconstructs the 
magnetization disappear as shown our divacancy results. 

Figures 5(b) and 5(c) shows the density of states for the 
(6,6)+2V CNT calculated with supercells of 12.46 and 
24.96 A length, respectively, which are compared with 
the density of states of the perfect (6,6) CNT [Fig. 5(a)]. 
We find that for the higher defect concentration the sys- 
tem has a metallic character with a small energy gap 
just above the Fermi energy, whereas for the lower con- 
centration it become a semiconductor with a energy gap 
of about 0.15 eV. This behavior is originated in a dis- 
persive CNT state which cross the Fermi level at about 
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FIG. 6: Density of states for the (10,0) CNT with a diva- 
cancy. (a) Perfect (10,0) CNT (a) (10,0)+2V calculated with 
a supercell of 12.96 A length, (b) (10,0)+2V calculated with a 
supercell of 25.92 A length. Dashed line indicate the position 
of the Fermi Level (Ef). 

0.8rX. In the more diluted defect concentration this state 
becomes flat opening a energy gap and turning the sys- 
tem a semiconductor. Similar results are found for the 
(10,0)+2V CNT shown in Fig. 6. Here, the perfect (10,0) 
CNT is a semiconductor with an energy gap of 0.8 eV 
[see Fig. 6(a)]. For the higher vacancy concentration 
[Fig. 6(b)] a defect state crosses the Fermi level at about 
O.irX while a CNT state rises close to the Fermi level, 
giving a metallic character to the system. However, for 
the more dilute defect concentration, the (10,0) +2 V CNT 
becomes semiconducting mainly due to the defect state 
above the Fermi level which becomes flat opening a en- 
ergy gap of about 0.15 eV as shown in Fig. 6(c). It is 
interesting to note that the (8,8) and (14,0) CNTs with 
a divacancy, calculated with supercells of L « 13 A, ex- 



hibit an energy gap of about 0.15 eV. This suggests that 
for more diluted defect concentrations in CNTs with di- 
ameters of at least 11 A, the gap energy should increase 
due to the flattening of the defect states close to the 
Fermi energy. 

V. SUMMARY AND CONCLUSIONS 

We have studied the monovacancy and divacancy in 
armchair and zigzag CNTs with diameters of about 8 and 
11 A and considering two concentrations of defects along 
the tube axis. Our results show that the CNTs with a 
monovacancy exhibits a ferromagnetic ordering induced 
by an undercoordinated C atom. The magnetic moment 
in the defect equilibrium geometries fluctuate among 0.3 
and 0.8 /ie, depending on the defect concentration and 
on the diameter of the CNTs. The electronic density of 
states along the CNTs axis show that the monovacancy 
does not change the electronic character of the tubes. On 
the other hand, once a divacancy is created the neighbor- 
ing C atoms spontaneously reconstruct around the defect 
forming an octagon with two adjacent pentagon in oppo- 
site directions. Because all the C atoms are three coordi- 
nated in this defect, CNTs with a divacancy are not mag- 
netic. However, the electronic properties of CNTs change 
when they contain a divacancy where both armchair and 
zigzag CNTs become semiconductor having a energy gap 
of about 0.15 eV. The above properties suggest the pos- 
sibility to achieve nanoscale magnets and magnetic pat- 
terns with potential applications in spintronics and bio- 
compatible magnetic materials. 
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